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ABSTRACT 

We constrain the distribution of calcium across the surface of the white dwarf star G29-38 by 
combining time series spectroscopy from Gemini-North with global time series photometry from the 
Whole Earth Telescope. G29-38 is actively accreting metals from a known debris disk. Since the 
metals sink significantly faster than they mix across the surface, any inhomogeneity in the accretion 
process will appear as an inhomogeneity of the metals on the surface of the star. We measure the flux 
amplitudes and the calcium equivalent width amplitudes for two large pulsations excited on G29-38 
in 2008. The ratio of these amplitudes best fits a model for polar accretion of calcium and rules out 
equatorial accretion. 
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1. INTRODUCTION 

The well known G29-38 is a bright white dwarf star 
with a H atmosphere (DA) that falls in the instabil- 
ity strip (DAV), has metals apparent in its photosphere 
(DAVZ), and has a known debris disk (DAVZd). The 
combination of these attributes affords us the opportu- 
nity to measure the geometry of accretion from the debris 
disk onto the white dwarf star. 

As a DAV, G29-38 shows multi-periodic, large ampli- 
tude pulsations (> 2%) with periods between 200 and 
1000 s. These are non-radial, g-mode pulsations that 
likel y arise from convective driving a.t the surface of the 
star (IBrickhilllll99"llGoldreich fc WdHoOflh . Each pulsa- 
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tion is a standing wave that can be described by three 
eigenvalues, fc, ^, and m. The first describes the ra- 
dial nodes while the last two, {I, m), denote the spher- 
ical harmonic, or the surface geometry of the pulsation. 
G29-38 does not always excite the same modes, select- 
ing several modes from the set of possible eigenvalues 
emman et all 119981 ). Matching the observed peri- 
ods with the correct eigenvalues is an important step 
in order to use pulsators like G29-38 to measure stel- 
lar interiors with a technique known as asteroseismol- 
ogy. As a bright DAV, G29-38 has been the target 
for developing techniques to measure {£, m) of each 
excited mode. Techniques such as time series spec- 
troscopv (Ivan Kcrkwiik ct al. 2000; Clemens ct al. 200^ 
iThompsonetTall 120081) . combination modes (fWul 12001). 
and light curve fitting (jMontgomervl [20051 ) have found 
excited £=2 modes among primarily £=\ modes. 

G29-38 became part of a larger mystery when the 
metals Ca, Mg and Fe were found in its photosphere 
(jKoester et al.lll997l ). White dwarf stars have very large 
gravities, and the diffusion timescale for Ca to precipi- 
tate out of the photospheres for those with temperatures 
less than ^ 20 , 000 K is much less than t heir cooling times 
(IKoestedl2009l:IKoester fc Wilkenll2006D . For G29-38 this 
diffusion timescale is less than 2 weeks, meaning G29-38 
must be actively accreting metals. The source of the 



debris disk ( 


Tokunaga et al.l 


iReach et al.l 


20051). The in- 



flfM IGraham et aT1fl99 
frared excess is well modeled by sub-micron silicates 
and amorphous carbon lying between 1-5 R0 from the 
star and is pos sibly formed from disrupted minor bod- 
ies (jFarihi et al...2009: Jura, 2008. '). The material accretes 
onto G29-38, producing the metal lines seen in its pho- 
tosphere. 

The accretion process onto G29-38 does not need to 
be geometrically uniform. For instance, the dust around 
G29-38 could be a thin disk that accretes material near 
the stellar equator. Alternatively, a magnetic field may 
guide material onto the star's poles. Once the metals 
are accreted they begin to sink out of the photosphere. 
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At the same time the metals are mixed across the stellar 
surface by the turbulent viscosity produced by the con- 
vection zone. As long as the timescale for gravitational 
settling is much shorter than the mixing timescale, the 
inhomogeneity of met als on the sur face of the star will 
be maintained. , Montgomery et all (2008) estimate the 
gravitational settling to be 30-50 times shorter than the 
horizontal mixing. If the accretion is indeed non-uniform, 
then the surface of G29-38 will have a non-uniform dis- 
tribution o f Ca acro ss its surface. 

Montgo merv et a l. (2008) proposed a technique to use 
the pulsations of G29-38 to measure any inhomogeneity 
of the surface metals. This technique is based on the fact 
that the strength of the metal lines is a function of tem- 
perature. As the pulsations cause temperature variations 
on the surface of the star, the equivalent width (EW) of 
any metal line will vary along with the pulsations. For 
G29-38, the size of thi s effect is as large as 10% for a 
single pulsation mode ()von Hippel fc Thompsonll2007l ). 
The pulsation temperature variations are described by a 
spherical harmonic, and since g-modes all have i>l, the 
temperature variations are not uniform across the sur- 
face. As a result, the amplitude of the EW variations 
relative to the photometric amplitude (/ = Aew/^Aux) 
will depend on the spherical harmonic of the pulsation 
mode and on the distribution of the metals. When the 
pulsation geometry aligns with the metal inhomogeneity, 
Aew will be larger. For example, if the metals followed 
an equatorial distribution, the EW variation due to the 
pulsations would be larger for a (1, 1) mode that pri- 
marily pulsates at the equator than a (1, 0) mode that 
primarily pulsates at the poles. 

The goal of this paper is to measure the geometry of Ca 
across G29-38 using the technique of iMontgomery et al.l 
((2008). We obtained time series spectroscopy from 
Gemini-North in order to measure the EW variations of 
the Ca K 3933A line. We also obtained time series pho- 
tometry from a collection of telescopes associated with 
the Whole Earth Telescope to accurately measure the 
photometric amplitudes of the pulsation modes. G29- 
38 excited two large amplitude triplets during our 2008 
data campaign, each split by 15 fiUz. We use the two 
largest modes, with different values of m, to constrain 
the inhomogeneity of Ca on the surface of G29-38. Fi- 
nally, we apply this same analysis to previous time se- 
ries spectroscopy of G29 -38: the VLT (Ve ry Large Tele- 
scope) data published bv lThompson et all ( 20081 and the 
Keck Observatory data published bv lvan Kerkwijk et al.l 
((2000h . 

2. DATA 

We received 9.5 hr of time series spectroscopy of G29- 
38 from GMOS (Gemini Multi-Object Spectrograph, 
iHook et al.l 120041 ) on Gemini-North in 20 08. Concur- 
rently we used the Whole Earth Telescope (jNather et al.l 
[l990fi to obtain over 300 hr of time series photometry 
from 10 different telescopes. The third priority ("band- 
3" ) , queue-scheduled Gemini observations were obtained 
during two nights in October and three nights in Novem- 
ber. These observations were sufficiently separated that 
we treated the data as two separate measurements of 
G29-38's pulsations. Luckily, the same large pulsation 
modes were present during both months, allowing us to 
easily average our measurements. 



2.1. Photometry 

We obtained the time series photometry for two dis- 
tinct purposes. First, we needed enough data to clearly 
identify the pulsation frequencies present in G29-38. Sec- 
ond, we needed to correctly measure the flux amplitude 
of each pulsation. We obtained photometry of G29-38 
from early August through early December in 2008. The 
majority of the coverage is in September, while the Gem- 
ini spectra were not taken until October and November. 
Several observers at small telescopes were on-call to take 
these observations when Gemini observed. The observa- 
tions we present here are split into three groups. The 
first, in September, is the most extensive but also lacks 
corresponding Gemini observations. The second, in Oc- 
tober, and the third, in November, each have ~5 hr of 
accompanying time series spectra. See Tables 1-3 for a 
list of the photometric observations. 

In September we obtained 29 separate runs from 8 dif- 
ferent telescopes, totaling more than 177 hr of photom- 
etry. In October we obtained 13 separate runs from 5 
different telescopes, totaling more than 50 hr of photom- 
etry. In November, we obtained 15 runs from 6 differ- 
ent telescopes, totaling more than 63 hr of photometry. 
Each observatory used 5-15 s exposure times with a mini- 
mal readout time of less than 15 s. We used a broadband, 
red cutoff filter, BG4O or similar, to reduce the variations 
between the sensitivities of each telescope and CCD cam- 
era. The PROMPT telescope (Panchromatic Robotic 
Optical Monitoring and Polarimetry Telescopes) , lacking 
a BG4O, used a Sloan g filter. Though the g filter lets 
through less red light, this has little effect on the mea- 
sured amplitudes of the DAV pulsation modes (<3%). 

To reduce the photometric images, we applied dark, 
bias and flat-field corrections to our image frames. 
We then performed aperture photometry using IRAF's 
PHOT routine on each image. We used the aperture 
that gave the highest signal-to-noise for each night of 
data. To create the final light curve, we divided by 
a bright star, removed bad points, and corrected for 
differential extinction using a program called WQED 
(|Thompson fc MuUallvl I2009D . Finally we applied a 
barycent ric correction to our times using the method of 
IStumpfB p9801 

2.2. Spectra 

We obtained a total of 9.5 hr of time-series spec- 
troscopy on 5 different nights in 2008 (October 16, Octo- 
ber 18, November 03, November 04, November 06), ap- 
proximately half of our allocated queue-scheduled Gem- 
ini observations. Each spectrum has an exposure time 
of 55 s and a readout time of 13 s, yielding a Nyquist 
period of 136 s. We reduce these spectra from Gemini 
IRAF GMOS routines, applying flats, biases and arc- 
lamps taken at the beginning of the night. The resulting 
dispersion is 0.9 A while the resolution is <3.2 A, set by 
a 1.5" slit. The signal-to- noise of the spectra is typi- 
cally 80, while the average spectrum for one night has a 
signal-to- noise of ~ 300 (see Figure[T]). The spectra cover 
the optical range from 3700-6400 A, including both Mg 
II (4480 A) and Ca K (3933 A). While both are present 
in the spectra, the Mg line is much weaker, and we are 
unable to measure variations in its line strength. 

We use the spectra primarily to measure the varia- 
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tions in the EW of the Ca hne. However the prominent 
Balmer hnes also change with the pulsations. We use the 
H lines to measure the line shape variations to identify 
the spherical degree, £, of each mode (see 

3. MEASURING THE PULSATIONS 
3.1. Flux Variations 

We perform a Fourier transform (FT) of the light 
curves from the September, October and November data 
sets (see Figure [2|). To judge the alias peaks created from 
the gaps in our data, we present the Fourier windows of 
each data set in Figure|3l G29-38 appears to have excited 
essentially the same modes during all our observations. 

We fit and remove the largest peaks found in t he FT 
using the program PERIOD04 (jLenz fc Breged l2005j ). 
While each group of data is dominated by the same sets 
of modes, their frequencies a ppear to wander from month 
to month. [Provencal et al.l (|2009f ) have observed similar 
behavior in other pulsating white dwarf stars. We use the 
September photometry as a guide to select the dominant 
modes and their combinations. We identify modes from 
subsequent months with the same name (e.g., Fl) if they 
lie within 3 /iHz. In the rare case that a statistically sig- 
nificant mode does not exist in the FT from that month, 
we perform a nonlinear fit at the largest peak near the 
September frequency (F4 in November is only marginally 
significant, but was fit for completeness). While signif- 
icant peaks may remain, we do not discuss them here. 
We use only the largest modes and their combinations to 
determine the metal distribution across G29-38. 

We provide the resulting list of frequencies, amplitudes 
and phases from each data set in Tables 4-6. The zero 
time for all phases is BJED=2454700.0. We calculate 
errors for the photometry with a 100 point Monte Carlo 
simulation using the program PERIOD04. We also fit 
the combination and harmonic peaks present in the FT, 
by constraining each frequency to the exact frequency 
derived by the sum of the parent mode(s). 

3.2. Equivalent Width Variations 

The convection zone of a DAV extends into its photo- 
sphere, so the time series analysis of its spectral lines is 
simpler than for other pulsating stars. The observed pul- 
sations are g-modes that are non-propagating ("evanes- 
cent" ) in the convection zone and the photosphere. Ad- 
ditionally, the short turn over time (< 1 s) of the convec- 
tion zone efficiently mixes the velocity variations, mak- 
ing them constant as a function of optical depth, and 
the total (i.e., bolometric) flux variations are also con- 
stan t through these outer layers (see Figures 2 through 
5 of iWu fc Go ldrcich 1999). From our models we esti- 
mate that these conditions are satisfied from the outer 
atmosphere down to optical depths of r > 1000. In ad- 
dition, convective mixing removes any depth dependence 
of the Ca abundance. Thus, any local region of the star's 
surface can be modeled with a static model atmosphere 
having the instantaneous local values of the Tcff, log 5, 
and Ca abundance. This is to be contrasted with the 
case of the roAp stars, in which the eigenfunctions can 
have sig nificant variations even at small optical depths 
e.g., see IMedupe fc IOart3[T99l . 

To measure the Ca EW variations, we normalize by 
the continuum and then fit a Gaussian to the Ca line in 



each spectrum. For low resolution data, simulations show 
us that fitting a Gaussian more accurately measures the 
area of the line than does numerical integration across 
the line. Since all the spectra were acquired with the 
same spectrograph and set-up, we averaged the spectra 
to determine the shape of the continuum. (The contin- 
uum in this case is actually the wings of the He and H8 
lines.) We fit a 5'*' order polynomial across 35 A of the 
spectrum, after removing the central 5 A of the line. This 
shape is used to flatten the continuum of each spectrum. 
Prior to measuring the EW, we normalize a second time 
by a linear trend to remove any remaining brightness and 
color variations due to the star or the Earth's atmosphere 
(again, we do not include the central 5 A of the line). By 
fitting the shape of the continuum with the average spec- 
trum we reduce the overall noise that would come from 
establishing the continuum with much noisier, individual 
spectra. The EW is computed as the area removed by 
the best fit Gaussian divided by the continuum level for 
each spectrum. 

To create fractional variations in EW, we divide by 
the average EW of the line from each night. We mea- 
sure the EW amplitudes by fitting these variations at 
the exact frequencies measured from the accompanying 
photometry. In Figure [3] we present the FT of the EW 
measurements. The noise level, measured from the aver- 
age power of the high frequency end of the FT, is 2.8% in 
the October data and 2.3% in the November data. Since 
we know the frequency from the photometry, the peaks 
found at the frequencies for Fl and F2 are significant. 
In both November and October, the probability that Fl 
in the EW FT is due to noise alone is 1/200, and simi- 
larly the prob ability that F 2 is due to noise alone is less 
than 1/2000 (|KepledlT99l . Tables 5 and 6 summarize 
the amplitudes of the EW variations along with the for- 
mal errors for the largest modes present in October and 
November. 

The above method for measuring the Ca EW suffers 
from systematic underestimates of the EW because of 
poor continuum normalization. In the next section we 
consider how the continuum fitting impacts our EW mea- 
surements and we quantify this effect. 

3.3. Systematics in Measuring the Continuum 

The EW amplitudes are measured relative to the av- 
erage EW. If we fit a continuum that mis-measures the 
average EW, then the EW amplitude itself will be wrong 
because we will have falsely added or subtracted area 
from the line. In general we find that we pick an average 
EW that is too small because the continuum includes the 
wings of the broad line, making the assumed continuum 
across the line have less flux than it should. 

We estimate this of fset by fitting model spectra con- 
taining the Ca K line (jWeidemann fc Koesterl[T98i ). We 
create >200 model spectra with a dispersion (0.9 A), 
resolution (3.2 A), and signal-to-noise ('^300), approx- 
i mately equivalent to our a verage Gemini spectrum 
(jWeidemann fc Koesterl 1198^ . After fitting the contin- 
uum and measuring the EW with the same procedure as 
the data, we find that the measured EW is only 70 ± 1% 
of the actual, model EW. If not corrected for this effect, 
the measured Ca EW amplitude would be too large. 

Our assumption that the continuum shape is identical 
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Figure 1. A portion of the individual and average spectrum of G29-38 taken with GMOS on Gemini-North. Each spectrum has been 
normalized with a low-order spline function for presentation purposes. The Ca K (3933 A) line has been labeled. 
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Figure 2. The Fourier transforms of light curves measured during the three G29-38 photometry campaigns. The frequencies of Fl, F2 
and their combinations are denoted at the top of the plot. 



throughout the pulsation also causes a sj'stcmatic offset, 
though a small one. Because we normalize all spectra 
by the same continuum, established from the mean spec- 
trum, there is a slight systematic linear trend to the off- 
set. Hotter spectra have weaker Ca lines than the cooler 
spectra. Additionally, a change in tempcratiirc alters 
the shape of the adjacent H lines and hence alters the 
effective continuum for the Ca line. Our technique un- 
derestimates the Ca EW in the hotter spectra to a larger 
extent than it does the cooler spectra. For example, the 
measured EW is 73% of the actual EW for measured 
EWs that are 10% larger than the mean. Based on this 
and other calculations, we correct our measurements us- 



ing the form 70% + 0.3a; where x is the per cent difference 
between the measured EW and the average EW. As the 
pulsations cause EW variations by at most 20%, this ad- 
ditional linear trend is small, but significant. 

Tables 5 and 6 give the EW amplitudes after correcting 
by these effects. We find the average Ca EW to be 279 ± 
10 mA. 

4. GEOMETRY OF PULSATION 

To compare our measiircmcnts of the ratio between the 
Ca EW and the flux amplitudes, /ca= ^EwMfiux (where 
we use the subscript to refer to either the element or the 
mode measured), we must first determine each pulsa- 
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Figure 3. The Fourier window for each month of data is shown 
at lOOO/^Hz. September had the most coverage and thus the fewest 
aUas peaks. 

tion mode's surface geometry. Simply put, we need the 
spherical harmonics, {£, m), and the inclination of the 
pulsation axis relative to our line-of-sight. Recently, sev- 
eral methods for determining the spherical degree and 
azimuthal order have been developed. Each method in- 
dividually leaves room for doubt, but used together we 
can reach an agreement on the geometry of the pulsation 
modes. We use the techniques of asteroseismology, com- 
bination modes, and line-shape variations to determine 
the {£, m) of the two largest pulsation modes. The incli- 
nation angle is more difficult to determine and is left as 
a free parameter when we com pare our measure ments to 
the inhomogeneity model of M ontgomerv "eTall (poos) . 

4.1. Multiplet Structure 

Stars pulsate at frequencies dictated by their interior 
structure. Asteroseismology, the technique of using the 
pulsation frequencies to determine the properties of the 
stellar interiors, predicts that when a star rotates, an £ — 
1 mode will be split into 2i+l modes with different values 
of m. The size o f the frequency spli tting depends on the 
rate of rotation (jWinget et al.lfTOQll) . Both of our modes 
are part of clear triplets split by 15 /iHz, giving us the 
likely guess that each mode is £=1. Figure [S] shows both 
triplets, labeled according to Table 4. Unfortunately, 
white dwarf pulsators do not always excite all members 
of a multiplet; an £=2 multiplet may also appear as a 
triplet leaving some doubt to the identification of our 
modes. If both are £ = 1 then Fl, as the first member 
of one triplet, and F2, as the central member of another, 
yield the most likely {£, m) identifications of Fl=(l, -1) 
and F2=(l, 0). 

4.2. Line shape variations 

Since the extent of limb darkening varies across the 
broad H features in white dwarf spectra, the shapes of 
the H lines pe riodically vary a c cordin g to the £ of the pul- 
sation mode. iRobinson et al.l ()1995f) first introduced the 
technique to measure how the amplitude of the pulsation 
varies at different wavelengths for pulsating white dwarf 
stars. Since then the technique has been successfully ap- 
pHed to bright DAV aiid DBV stars (e.g. iClemens et al.l 
[2000t IKotak "erall[200l [200l . The line shape variations 
are independent of the inclination angle and m, a nd can 
only be used to measure £. Recently, iThompson et al.l 



(|2004 12008D have shown that fitting a pre-specified func- 
tional form to the spectral lines allows this technique to 
be applied to poorer q uality data. Since our ob servations 
are of the same star as IThompson et aD ([2008!), with ex- 
actly the same technique we can directly compare our 

results to the previous time series spec troscopy. 

Foll owing the technique described by IThompson et"al] 
(12001 . we measure spherical degree by fitting the H/3 
and H7 lines in our Gemini spectra. We describe the 
variation in the shape of the Balmer lines by fitting the 
combination of a Gaussian and a Lorentzian to each spec- 
tral line. We measure the EW of both the Gaussian and 
the Lorentzian and measure their variations at the known 
pulsation periods to determine the fractional amplitudes 
(denoted Gevv and Lewj respectively). We then normal- 
ize those amplitudes by the fiux amplitude of the mode 
measured by the fitted continuum at the center of the 
Balmer line (e.g., ^4862)- By plotting the amplitude of 
the variations in the normalized areas of the Lorentzian 
against the Gaussian, we can distinguish between differ- 
ent values of £. While there is a known offset between 
the models and the observed modes, different spherical 
degrees appear to follow the basic trend outlined by the 
models. 

We limit this analysis to the two largest photometric 
modes (Fl and F2) as they are the only modes large 
enough to show clear line shape variations. Figure [6] 
shows the H7 and H/3 variations compared to both a 
model appropriat e for G29-38 and the previous observa- 
tions of this star (jGlemens et al.l[2000l IThompson et al.l 
l2008f ). Our Gemini observations were of high enough 
quality to measure Fl in October and both Fl and F2 in 
November. The formal error bars from fitting the Gaus- 
sian and Lorentzian functions at the known periods are 
represented in Figure [6l Both Fl and F2 lie amongst 
the previously observed £=\ modes, making it the most 
likely interpretation of the data. However, given the size 
of the error bars and the ambiguity in where an £—1 
mode should lie, there is still a nonzero chance that Fl 
is £=A and F2 is £=2. 

4.3. Combination modes 

We use combination modes to identify (£ , m) of th e 
pulsation modes in the manner described by |^ (|200l . 
The ratio of the combination amplitude ( Ajj ) to the par- 
ent amplitudes (Aj and Aj) is given by Rr=AyJn\\ A\A\, 
where riij is 1 for harmonics and 2 for combinations. iWul 
using the surface geometry of each mode and the 
inclination angle, makes a prediction for this value. Her 
analysis only uses the first order perturbation, and thus 
the mode identification technique is most successful when 
applied to lo w-amplitude pulsators (Handler et al. 2003 
lYeates et aLll2005( ). While our Fl and F2 are relatively 
large amplitude modes, we find that this combination 
analysis yields consistent results and is still useful for 
G29-38. 

Figure [7] shows the predicted Rc values along with 
those values measured for the combinations of Fl and 
F2 on G29-38. We average the values over the three 
photometric runs (September, October and November) 
and use a standard deviation of the three measurements 
to represent the error bars (the formal error bars are sig- 
nificantly smaller). The harmonic frequencies are most 
valuable in mode identification because the two parent 
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Figure 4. The Fourier transforms of tiie fractional Ca EW variations during October and November. Fl and F2 are labeled at the top 
for reference. 



modes have the same {£, m). The harmonic 2F2 hes sig- 
nificantly highe r than the harmonic 2F1. This is only 
consistent with IWul (|2001D when F2 is an (1, 0) mode 
and Fl is (1, -1). The combination theory then predicts 
the amplitude of a (1, -1) combined with a (1, 0) to be 
the same amplitude as the harmonic of (1, F1-I-F2 
thereby supports the results of the harmonics. This the- 
ory clearly states Fl=:(l, -1) and F2=(l, 0), as suggested 
by the triplet structure. 

4.4. Oblique Pulsator Model 

In certain types of pulsating stars, such as /? Cephei 
and roAp stars, the magnetic field can cause the pul- 
sations t o be aligned with an axis other than the rota- 
tion axis (iShibahashi fc Aert£|[2000t iKurtz fc Shibahashil 
Il986t lKurtdll982f ). In this oblique pulsator mode, a sin- 
gle pulsation will be split into 2.^-1-1 frequencies. These 
extra frequencies do not represent pulsation modes of the 
star, instead they are caused by a pulsation, not aligned 
with the rotation axis, rotating around the star. If this 
is the cause for the triplets in our data of G29-38, the 
15 /iHz splittings of our two modes indicate a rotation 
p eriod of ^18.5 hr. 

[Montgomery et al] (|2010l ) offer the first evidence of 
oblique pulsation on a pulsating white dwarf star and 
so we mention this possibility for completeness. If this 
model applies to G2 9-38 , the asteroseismological argu- 
ments presented in i|4.1l are not applicable. However, 
projecting the off-axis pulsation onto the rotation axis 
would yield triplets with exactly the same (£, m) ge- 
ometry as determined from seismology, except that now 
the frequency splittings are given by the inverse of the 
rotation period. Thus, we arrive at the same mode iden- 
tifications in either case. 

^® The 1=2 harmonics lie significantly off the plot in Figure [7] 
for any inclination angle, except for (2, 1), which lies near Rc ~ 3. 
However if Fl is (2, 1), the amplitude of F14-F2 should be much 
larger than observed. 



For a detailed analysis of the Ca geometry incorporat- 
ing this oblique pulsator model, we would need to make 
refinem ents to the model presented bv lMontgomerv et al.l 
(|2008f l. However, our data are not sufficient to make any 
but the most general conclusions, so we leave this exer- 
cise to future papers should the oblique pulsator model 
be deemed appropriate for G29-38. 

5. MAPPING CALCIUM ACROSS G29-38'S SURFACE 

We compare our mea surements of /r;a.= ^EwM fiiiy to 
the model presented bv lMontgomerv et al.l (|2008l ) for Fl 
and F2; see Tables 5 and 6 for the Jca values for Fl and 
F2 in October and November. Because we see the same 
modes in both October and November, we take the av- 
erage of the two measurements. Fl, as a (1, -1) mode, 
has an average /ca value of 1.7 ±0.3, and F2, as a (1, 0) 
mode, has an average /ca value of 2.9 ± 0.5. Figure [S] 
presen ts a direct comparis o n of t hese two measurements 
to the iMontgomerv et aH (|2008f ) models (we note that 
the predicted /ca values for (1, 1) are identical to those 
predicted for (1, -1)). Next, we discuss which Ca geom- 
etry is most likely and then consider previous data sets 
from the VLT and Keck observatories. 

5.1. Comparison to the models 

First we consider the simplest case, a uniform accretion 
model. Here the different geometries of the m eigenval- 
ues are inconsequential because the Ca has a uniform ge- 
ometry. A ll measurements of Ica. fo r the same £ should 
be equal. IMontgomerv et all ()2008f ) calculate that /ca 
should be 2.7 for t=\. Comparing our two measurements 
to this model yields a Xr = 4.3, or less than a 1% chance 
that this is the appropriate model. Because of our inabil- 
ity to accurately measure the continuum when measuring 
the EW, these two values may both suffer from the same 
fractional offset in EW amplitude and thus /ca- Even 
with the presence of such a systematic effect, the ratio 
of these two measurements would be 1 if the uniform 
model is appropriate. Instead, the ratio of their /ca val- 
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Figure 5. A schematic illustration of the two 15 /^Hz triplets observed in 2009 on G29-38, each panel represents a different month of data. 
Each arrow represents the frequency and the amplitude of a nonlinear fit to the complete lightcurve of each month. The FT of the residuals 
is shown in black for each month. 
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from previous G29-38 spectra are represented as circles Ijvan Kerkwiik eTan[2000l : [Thompson et al.|[200a i. The previously identified ^=2 
modes are presented as x's. 



ues is /fi//f2 = 0.6 ± 0.1, indicating less than a 1% 
chance these measurements agree. Primarily, the uni- 
form model fails to match Fl though it agrees quite well 
with F2. With only one outlier from two measurements, 
our data cannot completely rule out the uniform model, 
though they suggest the surface is inhomogeneous. 

Both the equatorial and polar models are more com- 
plex because the /ca values are dependent on the extent 
of the inhomogeneity, the inclination angle, and {£, m). 
For Fl and F2, three separate techniques agree with the 
conclusion that Fl is (1, -1) and F2 is (1, 0). This leaves 



only the inclination angle an d the extent of the inhomo- 
geneity. iMontgomervl (|2005f ) has measured the inclina- 
tion angle to be 65° by fitting the light curve of G29-38 
and we find an inclination of 55° in §4.3l bv using the am- 
plitude of the harmonic of F2. However, we do not rely 
on these measurements for our results and only use the 
inclination measurements to suggest that the inclination 
angle is not at the extremes. 

We co mpare our observations to the model described 
by M ontgomerv et al. (2008). where the rate of gravita- 
tional settling greatly exceeds the horizontal mixing rate 
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Figure 7. The ratio between the combination mode and parent 
mode amplitudes (Rc) is plotted against the frequency of the com- 
binations. Here we consider the harmonics of Fl and F2 along 
with their com bination. The predicted values of Rc (Wu, 200li; 
lYeates et al.l2005i ) for £=1 modes are shown (blue as the (1, 0) har- 
monic (at inclinations of 45°, 55° and 65°) and red as the (1, ±1) 
harmonic). The combination Rc values represent the average of 
the three photometric runs discussed in this paper, the error bar 
represents the scatter in those three measurements. The £=2 com- 
bination modes have large Rc values, except for the harmonic of 
(2, 1) which lies neariJc=3. This linear combination theory is con- 
sistent with the identification of F2 as (1, 0) and Fl as (1, -1). 

by turbulent viscosity. However if the horizontal mixing 
rate were increased, or the location of the accretion is 
wider than assumed, the extent of the inhomogeneity for 
either the polar or equatorial model would decrease. For 
the equatorial model we show in Figure [SI the ratio of 
these values, /(i,i)//(i,o): is between 2.3 and 3.7 depend- 
ing on the inclination of the pulsation axis. The polar 
case is the exact opposite; ^(1,0) '^iU always be larger than 
1(11), with the ratio lying between 0.22 and 0.38. In both 
cases the ratio will approach one if the inhomogeneity is 
not as extreme as our model assumes. 

Figure H] shows that the polar model fits better than 
the equatorial model because /fi//f2 = 0.6 ± 0.1 < 1. 
The Ica values for Fl and F2 differ from the equatorial 
accretion model by >lla. However, they differ from the 
polar accretion model by 1.5-2.5cr, depending on the in- 
clination angle. In both cases we could improve the fit 
by adjusting the amount of the inhomogeneity, i.e., cre- 
ating a larger region of Ca. Only a minor adjustment is 
needed to find a polar model that agrees with our mea- 
surements. There is no way to adjust the models to find 
agreement with the equatorial model. 

While we cannot entirely eliminate the uniform distri- 
bution, based on the two measurements of Fl and F2, 
we conclude that the Ca on G29-38 has a higher concen- 
tration at the poles than it does at the equator. 

5.2. Previous observations 

Keck and VLT time series spectroscopy also exist for 
G29-38, primarily taken for the purpose of i i denti- 
fication (Ivan Kerkwiik et all 120001: Iciemens et"all 120001: 
[Thompson et al.ll2008D ." For both, no adjoining photom- 
etry was taken, so all measurements of the pulsation flux 
amplitude are from the spectroscopy alone. While the 



photometry from the spectra is quite good, the frequency 
resolution is poor. The Keck and VLT spectral series are 
only 4 and 6hr in length, respectively. This does not 
allow us to distinguish splittings as close as the 15 /iHz 
splittings that we see in the Gemini observations. As a 
result we cannot assume that our measurements are not 
influenced by unresolved modes. Keeping this in mind, 
we measure the Ica values of the largest modes in the 
1996 Keck data and the 1999 VLT data in the manner 
explained above. 

We use t he reduced data s ets o f [van Kerkwiik et all 
(|2000D and [Thompson et al.l ([2008[ ) to measure the Ca 
EWs. The pulsation modes and the photometric ampli- 
tudes are taken directly from the respective papers. We 
measure the Ca EW for each individual spectrum as de- 
scribed in H3.2I As with the Gemini data, we likely under- 
estimate the continuum, causing overly large values of Ca 
EW amplitudes. Using model white dwarf spectra emu- 
lating the VLT and Keck observations, we calculate the 
appropriate corrections as described in We can only 
measure the Jca values for the largest two modes in each 
data set, providing us four more measurements of the in- 
homogeneity of Ca on G29-38. Our average Ca EW for 
the Keck and VLT data after accounting for these off- 
sets is 198 ± 10 mA and 298 ± 10 mA, respectively. The 
Keck data have a signiflcantly lower Ca EW than the 
VLT or Gemini observations, a fact previously noted by 
[von Hippel fcT hompson (2007). 

For these previous data sets to yield Ca inhomogeneity 
we need to know the {i, m) of the four modes. The line 
sha pe variations n i easure d by [Thompson et al.[ (|2008[ ) 
and [Clemens et al.l (|2000( ) tell us that all 4 modes are 
1=1. The values for m are less constrained. We can use 
the combination frequencies, as we did in ij4.3l to give us 
good guesses for m. Both the Keck and VLT data have a 
615 s mode. The amplitude of the harmonic for the 615 s 
mode is low (i?c~ 4) in both the Keck and VLT data; by 
looking at Figure[7l this indicates that it is a (1, 1) mode. 
lMontgomerv[ ([2005f ) also found this mode to be (1, 1) by 
fltting G29-38's light curve. The other large amplitude 
Keck mode is at 817 s, and the other VLT mode is at 
810 s. In both cases the amplitude of the harmonic is 
large {Rc ^ 13), like an (1, 0) harmonic should be. We 
tentatively conclude these modes are (1, 0). See Table[7| 
for our EW and /ca measurements along with the best 
guess at (^, to) for the 4 modes we measured in the Keck 
and VLT data. 

Figures [9| and [TOj show a comparison of the VLT and 
Keck data to the polar , equatorial and uni f orm ac cretion 
models calculated by [Montgomery et aP (|2008[ ). The 
VLT values for /ca are very similar to the Gemini obser- 
vations. The (1,1) mode is signiflcantly smaller than the 
(1, 0) mode, ^(i,i)/^(i,o) = 0-5 ± 0.2. Again, this cannot 
rule-out the uniform distribution, but it disagrees with 
the equatorial model. The two Keck modes, however, 
have identical values for Jca- If these modes do indeed 
have different m values, then the Keck data agrees best 
with a uniform distributiorF^. 

Another indication that the 1996 inhomogeneity is dif- 



IM ontgomerv et al.' f200J^) used measurements from 
Ivon Hipp el &; Thompson (2007) that had not been corrected 
for offsets introduced by mis-measuring the continu um a nd thus 
had a much higher value for the 615s mode (see ^3.3p . 
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Figure 8. The measured /ca values from the Gemini observations compared to the calculations of Mo ntgomerv et al.l lj200SI 1 as a function 
of inclination angle (i). The plot on the left depicts the models for the (1, 0) mode while the plot on the right depicts the (1, 1) mode. The 
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Figure 9. The measured value for /ca using the 1999 
VLT observations compa red to the model values presented by 
IMontgomerv et al.l II2008I ) as a function of inclination. The plot 
on the left depicts the models for the (1, 0) mode while the plot 
on the right depicts the (1, 1) mode. The grey box represents the 
1(7 error bars for the 615 s and 810 s modes. The measurements 
disagree with the equatorial model, see i|5.2l 

ferent than the 1999 epoch is the 2(7 difference in the 
/ca measurement for the 615 s, (1, 1) mode. A change 
in /ca indicates that the location of the Ca on the sur- 
face of G29-38 changed between 1996 and 1999. Except 
for the assumption that both data sets contain the same 
mode, this conclusion that the Ca geometry changed is 
independent of any model. 

5.3. Inhomogeneity from Velocities 

The g-mode pulsations cause the atmosphere to flow 
periodically across the surface of the star, resulting in a 
periodic Doppler shift of the Balmer lines at the period 
of the stellar pulsations. Due to the rapid mixing by the 
turbulent motions which extend into the photosphere, 
the pulsation velocities are independent of optical depth. 
Also, the amplitudes of the velocities depend on the ge- 
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Figure 10. The measured value for /ca using the 1996 
Keck observations compared to the model values presented by 
Montgomer v et al.l 1 12008 ) as a function of inclination. The plot 
on the left depicts the models for the (1, 0) mode while the plot 
on the right depicts the (1, 1) mode. The grey box represents 
the 1(T error bars for the 615 s and 817 s modes on the (1, 1) and 
(1, 0) plots respectively^ The data most closely match a uniform 
distribution of Ca, see ^5.21 



ometry of pulsa tion, so we can use them t o find different 
values of £ (see Ivan Kerkwiik et "all I2000D . In addition, 
the geometry of the velocity fields can indicate if any 
surface inhomogeneity exists. For instance, if metals lie 
inhomogeneously across the surface of G29-38, the ampli- 
tude of their pulsational velocities would not necessarily 
equal those measured from the H lines. The geometry of 
the velocity variations for these g-modes is proportional 
to the derivative of the flux variations. Thus, a larger ve- 
locity amplitude from the metal line compared to the H 
lines could result if the metals were concentrated where 
the velocities are the largest, near the flux nodal line, 
e.g., the equator of an (1, 0) mode. 

We fit the H/3 and H7 lines as described in 
[Thompson et al.l (|2004) . with a combination of a 
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Lorentzian and a Gaussian to measure the periodic 
Doppler shifts in the Gemini spectra. The Ca K hne 
is much smaller and we only fit a Gaussian to measure 
its velocity. We then create velocity curves for both the 
average H and the Ca K lines. Similarly, we measure 
the Ca K lines for both the Keck and VLT data sets; 
their H line amplitudes h ave a lready been publi s hed b y 
Ivan Kerkwiik et all (|2000[ ) and [Thompson et al.l (|2008[ ). 

We measure the velocity at Fl and F2 for the H ve- 
locity curves from the Gemini data (see Tables 5 and 6). 
For the H line, we also measure the ratio between the ve- 
locity and flux amplitudes (i?v — (^v/27r/)/Afiux)- Our 
velocity amplitudes agree with previous measurements of 
i=\ modes on G29-38, further evidence that Fl and F2 
are t=\ for our modes. We were unfortunately unable to 
measure the velocity amplitudes of the Ca K line from 
the 2008 Gemini and the 1999 VLT data. Both show no 
significant peaks in the FT with an uninteresting limit of 
llkms^^. This is significantly larger than the H veloc- 
ity amplitude and thus we can conclude nothing about 
the Ca inhomogeneity for these two data sets from the 
velocities alone. 

However, the 1996 Keck Ca velocity FT shows a 
significant peak for the 615 s mode with a velocity 
of 24±6kms~^. This measured velocity amplitude is 
4.4±1.3 times larger than the H velocity amplitude 
(5.4±0.8kms"i) measured bv lvan Kerkwiik et al.l (|20000 
for the same mode. It is also worth noting that another 
peak is present in the Ca velocity FT, at 7505 /iHz with 
an amplitude of 22±6kms~^. This 133 s mode does not 
appear in the Keck photometric FT and may indicate 
that these Ca velocity measurements are actually domi- 
nated by an unknown systematic effect. If the 615 s Ca 
velocity peak is real, the larger Ca velocity amplitude 
indicates that Ca is located on the surface in a way that 
maximizes the line-of-sight velocities. 

Precise modeling of the velocities that include the in- 
clination angle is needed to see how the Ca could be dis- 
tributed in order for its velocity amplitude to be much 
larger than the H velocity amplitude. As an example of 
how this might occur, the 615 s mode is (1, 1) and the 
largest line-of-sight velocities will come from the limb of 
the star for high inclinations. If all the visible Ca were 
located in the limb of the star during the Keck obser- 
vations, only the largest, line-of-sight velocities would 
be observed for Ca, while the H line would include re- 
gions with no line-of-sight velocity shifts. At the very 
least, this difference in velocity amplitude is an indica- 
tion that Ca is inhomogeneous in 1996, when the Keck 
observations were obtained. However the difference in 
velocities contradicts the /ca measurements for the Keck 
data which favor a uniform distribution of Ca. Possibly 
more complicated scenarios including spots of Ca could 
reconcile these two measurements. 

6. DISCUSSION OF ACCRETION 

If we consider our results for each epoch of measure- 
ment (1996, 1999 and 2008), then something about the 
accretion has changed in 1996 (Keck) when compared 
to 1999 (VLT) and 2008 (Gemini). The 1999 and 2008 
results indicate that the majority of the Ca accretes at 
the poles. This is most likely a signpost of magnetically- 
infiuenced accretion. However, in 1996 the accretion ap- 
pears to be different; the /ca values indicate a uniform ac- 



cretion. Perhaps the unusual re sults from the 1996 epoch 
should not be a surprise as Ivon Hippel fc ThompsonI 
(|2007D already noted a change in accretion at this time. 
The Ca EW was 70% smaller than in 1999 and the Ca 
EW showed fiuctuatio ns on tim e scales as short as 2 
weeks (but also see iDebes fc L6pez-Morales,200 8). Once 
again, the Keck data show that the accretion of G29-38 
was not entirely constant. 

These variations could be used to better understand 
either the disk or the accretion process. Since the accre- 
tion geometry indicates that the magnetic field was not 
infiuencing accretion in 1996, exactly when the EW was 
smallest, any Ca EW modulation may be connected to a 
modulation in the magnetic field. 

We do not know the magnetic field strength at the sur- 
face of G29-38. The non-detection of Zeeman splitting 
in this star only lim its the magnetic field to be < 10 kC 
(jKoester et al.lll998[ ). Yet the apparent detection of po- 
lar accretion implies magnetic fields. If G29-38's mag- 
netic field is strong enough and has a complex, evolving 
structure, then this magnetic field could modulate the 
quantity of gas that is accreted onto the star. In this 
case, the variations in surface Ca abundance would in- 
dicate that the timescale for significant changes in the 
magnetic field structure is less than two weeks. The 
plausibility of this scenario is, at present, limited by our 
knowledge of the field strength, how it infiuences accre- 
tion rates, and the physical mechanisms that would cause 
G29-38's magnetic field to change. 

This time scale of two weeks places constraints on the 
magnitude and global geometry of the proposed field. If 
the field has a large scale structure (say dipolar) and is 
non-zero outside of the star, then the energy contained 
in such a field is approximately 

E « {B^)V, (1) 

where V is the volume of the star and (i?^) is the value 
of averaged over the stellar surface. For a 10 kC field 
this is ^ 10^^ ergs. This is much larger than the total 
energy radiated into space by the star in a two week 
period, implying that building such a field on this time 
scale is not possible. On the other hand, if the field is 
less global and is non-zero only in the surface layers, then 
the energy in a 10 kG field would require only about 
1% of the energy radiated by the star in a two week 
period, which is energetically feasible. Thus, we may be 
seeing the effect of smaller-scale, local variations in the 
magnetic field which nevertheless are able to perturb the 
accretion geometry. 

7. CONCLUSIONS 

We combined our Gemini spectroscopy and Whole 
Earth Telescope photometry to constrain the accretion 
geometry of Ca on G29-38. First, we identified the pul- 
sation modes using a variety of techniques that took ad- 
vantage of all the information contained in the time-series 
spectroscopy and photometry. Second, we measured the 
Ca EW pulsation amplitudes and compared them to the 
flux amplitudes. This combination of measurements al- 
lowed us to constrain the geometry of C a on the surface 
of G29-38. These data from 2008 somewhat suggest the 
Ca is inhomogeneous: a polar distribution of Ca is most 
likely and the equatorial distribution has been ruled out. 
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If the Ca does lie at or near the poles, the accretion 
must be occurring there. This is the first evidence that 
a magnetic field is influencing the accretion process onto 
G29-38. 

We have shown how this technique may be used to 
directly measure the distribution of Ca on the surface of 
a white dwarf and thereby infer the accretion geometry; 
however a few improvements would have made our results 
more conclusive. Since we were allocated band-3 Gemini 
time, we chose a lower resolution grating to increase our 
odds of getting Gemini observations. A higher resolution 
would have improved our Gaussian fits of the Ca line and 
also helped establish the continuum across the line. Our 
photometry was significantly better in September, and 
had our spectroscopy come at the same time, we would 
have reduced our flux amplitude error bars by 30%. The 
best way to improve the strength of our conclusions is 
to simply take more measurements. With approximately 
twice as miich time series spectroscopy we woiild have 
reduced the error bars by 40% and then we could have 
convincingly distinguished between a polar and a uniform 
model. 

While previous time-series spectroscopy on G29-38 did 

not have the frequency resolution of our data, these data 
both confirm our polar conclusions (VLT observations) 
and present the tantalizing possibility that the accretion 
geometry may change (Keck observations) over time. As 
the Keck observations also coincide with a lower appar- 
ent Ca abundance, we postulate how the accretion geom- 
etry and rate may be related. With further observations 
of G29-38 we may establish that accretion in this white 
dwarf debris disk system is not only influenced by mag- 
netic flelds but is also dynamic in nature. 
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Table 1 

Photometric observations of G29-38 in September 2008. 



Observatory 


Date 






Start (UTC) 


Length (hr) 


OPD Brazil 


2008- 


■09- 


■20 


03: 


;55; 


;17 


3; 


;12; 


;21 


Apache Point Im 


2008- 


■09- 


■21 


05; 


;53; 


;39 


5; 


;05; 


;41 


Meyer Observatory 0.6m 


2008- 


■09- 


■21 


06; 


;17; 


;08 


4; 


;14 


;50 


Meyer Observatory 0.6m 


2008- 


■09- 


■22 


03: 


;01: 


;27 


6: 


;59: 


;50 


OPD Brazil 


2008- 


■09- 


■23 


00: 


;37: 


;13 


5: 


;32; 


;29 


McDonald Observatory 2.1m 


2008- 


■09- 


■23 


02; 


;11; 


;03 


9; 


:27; 


;25 


Meyer Observatory 0.6m 


2008- 


■09- 


■23 


02: 


;23: 


;40 


4: 


;40: 


;39 


Apache Point Im 


2008- 


■09- 


■23 


04; 


:50; 


;44 


6; 


;04; 


;59 


OPD Brazil 


2008- 


■09- 


■24 


00: 


;18: 


;02 


6: 


;56: 


;20 


McDonald Observatory 2.1m 


2008- 


■09- 


■24 


03; 


;10; 


;58 


2; 


:00; 


;25 


Meyer Observatory 0.6m 


2008- 


■09- 


■24 


03; 


;13; 


;03 


7; 


;29; 


;50 


Apache Point Im 


2008- 


■09- 


■24 


07; 


;47; 


;40 


2; 


;28; 


;42 


McDonald Observatory 2.1m 


2008- 


■09- 


■24 


08; 


;13; 


;13 


2; 


;51; 


;20 


McDonald Observatory 2.1m 


2008- 


■09- 


■25 


02; 


;06; 


;36 


6; 


;07; 


;00 


Meyer Observatory 0.6m 


2008- 


■09- 


■25 


04; 


;13; 


;00 


5; 


;59; 


;50 


Mcdonald Observatory 2.1m 


2008- 


■09- 


■26 


02; 


;04; 


;08 


8; 


;14 


;15 


Meyer Observatory 0.6m 


2008- 


■09- 


■26 


06: 


;27: 


;05 


3: 


;59: 


;49 


Apache Point Im 


2008- 


■09- 


■26 


07; 


;42; 


;35 


3; 


:03; 


;31 


Tuebingeii Observatory 


2008- 


■09- 


■26 


19: 


;12: 


;43 


7: 


;18: 


;46 


McDonald Observatory 2.1m 


2008- 


■09- 


■27 


02; 


;02; 


;32 


7; 


;22; 


;05 


Meyer Observatory 0.6m 


2008- 


■09- 


■27 


02: 


;02: 


;55 


7; 


;59: 


;50 


Mt. Suhora 0.6m 


2008- 


■09- 


■27 


18: 


;37: 


;59 


5; 


;35; 


;00 


Krakow 0.5m 


2008- 


■09- 


■27 


18: 


;53: 


;04 


2; 


;17; 


;10 


Tuebingen Observatory 


2008- 


■09- 


■27 


20; 


;30; 


:24 


6: 


;25: 


;36 


McDonald Observatory 2.1m 


2008- 


■09- 


■28 


02; 


;06; 


;17 


9; 


;07; 


;10 


KraJiow 0.5m 


2008- 


■09- 


■28 


17; 


;33; 


:49 


7; 


:06; 


;41 


McDonald Observatory 2.1m 


2008- 


■09- 


■29 


02; 


;57; 


;47 


7; 


;51; 


;55 


Meyer Observatory 0.6m 


2008- 


■09- 


■29 


07; 


:31; 


:01 


2; 


:59; 


;50 


Meyer Observatory 0.6m 


2008- 


■09- 


■30 


04; 


;19; 


:31 


4; 


:59; 


;50 



Table 2 

Photometric observations of G29-38 in October 2008. 



Observatory 


Date 






Start(UTC) 


Length (hr) 


Krakow 0.5m 


2008- 


■10- 


■05 


21:11:58 


2:10:43 


Mt. Cuba 0.6m 


2008- 


■10- 


■08 


01:51:51 


4:05:58 


Krakow 0.5m 


2008- 


10- 


■08 


18:35:56 


3:23:58 


Krakow 0.5m 


2008- 


■10- 


■11 


18:00:54 


5:44:32 


Apache Point Im 


2008- 


■10- 


■17 


01:23:51 


7:50:31 


Mt. Cuba 0.6m 


2008- 


■10- 


■17 


03:50:02 


2:50:14 


Mt. Suhora 0.6m 


2008- 


■10- 


■17 


17:29:49 


1:18:12 


Apache Point Im 


2008- 


■10- 


■18 


02:36:09 


6:32:30 


Gemini North 8m 


2008- 


■10- 


■18 


06:25:34 


2:17:02 


Krakow 0.5m 


2008- 


■10- 


■18 


16:45:29 


7:06:00 


Krakow 0.5m 


2008- 


■10- 


■19 


18:36:40 


1:06:00 


Mt. Cuba 0.6m 


2008- 


■10- 


■20 


00:03:13 


3:58:43 


Mt. Cuba 0.6m 


2008- 


■10- 


■21 


01:41:57 


4:40:41 
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Table 3 

Photometric observations of G29-38 in November 2008. 



Obscrva,tory 


Date 






Start(UTC) 


Length (hr) 


OVJ-rt.-t\, 4111 


zUUo- 


1 1 
■il- 


no 


Uv5.00.UU 


l.Oi.Uo 




onnc 
zUUo- 


■il- 


-Uo 


no. OQ . on 


0.17. Kn 
0. i ( .ou 


Vienna O.Sni 




1 1 
■ii- 


-Do 


ID.OU. ii 


D.ZO.ZO 


Mt. bunora U.om 


2008- 


■11- 


-03 


lo:o ( :zU 


5:33:02 


Apache Point Im 


ZuVo- 


1 1 
■il- 


r\A 

-U4 


ni .oo.oQ 
ui . zy .Zo 


D.Uz.Oy 


PROAAPT n Am 


zuuo- 


■11- 


DA 


uz . zo . zy 


o.Qft.i n 

Z.OU. J-U 


Mt. Suhora 0.6m 


2008- 


■11- 


-04 


17:19:39 


2:40:02 


PROMPT 0.4m 


2008- 


■11- 


-05 


01:37:41 


4:13:26 


Mt. Suhora 0.6m 


2008- 


■11- 


-05 


17:18:47 


6:02:12 


Apache Point Im 


2008- 


■11- 


-06 


01:55:29 


5:38:22 


Apache Point Im 


2008- 


■11- 


-07 


01:35:45 


5:27:49 


Apache Point Im 


2008- 


■11- 


-11 


02:53:48 


4:37:43 


Krakow 0.5m 


2008- 


■11- 


-11 


16:50:38 


5:44:40 


Apache Point Im 


2008- 


■11- 


-16 


03:10:10 


3:47:56 


Apache Point Im 


2008- 


■11- 


-17 


00:58:56 


5:53:54 



Table 4 

Measured variations in brightness of G29-38 from September 2008. 



mode 


frequency 


Aflux 


To, flux 


Rc 






% 


s 




Fl 


1209.610±0.006 


3.04±0.02 


742±1 




F2 


1112.443±0.007 


2.15±0.03 


779±2 




F3 


1224.701±0.009 


1.91±0.02 


579±2 




F4 


1239.721±0.014 


1.32±0.03 


654±2 




F5 


1127.520±0.016 


1.17±0.02 


142±4 




F6 


1097.446±0.017 


1.03±0.03 


791±4 




F7 


1293.038±0.023 


0.78±0.02 


561±4 




F8 


1522.498±0.019 


0.93±0.03 


144±3 




F9 


2502.323±0.013 


1.31±0.03 


291±1 




FIO 


1999.767±0.028 


0.63±0.02 


154±3 




Fll 


1015.214±0.031 


0.59±0.36 


890±291 




=2F1 


2419.220±0.000 


0.41±0.03 


412±4 


4.4±0.3 


=2F2 


2224.886±0.000 


0.76±0.02 


38±2 


16.3±0.6 


=F1+F2 


2322.053±0.000 


0.35±0.02 


422±5 


2.7±0.1 


=F1+F3 


2434.311±0.000 


0.24±0.02 


350±7 


2.1±0.1 


=F2+F3 


2337.144±0.000 


0.50±0.03 


383±4 


6.1±0.2 


=F2+F4 


2352.164±0.000 


0.34±0.02 


368±5 


6.0±0.2 


=2F4 


2479.442±0.000 


0.13±0.02 


75±16 


7.4±1.5 


=F3+F4 


2464.422±0.000 


0.16±0.02 


253±11 


3.1±0.2 


=F1-|-F6 


2307.056±0.000 


0.19±0.03 


9±9 


3.0±0.2 



Note: Variations are measured using the form Asin(27r/(t + To)). All To values are measured from BJED=2454700. Rc is the ratio of the 

parent to harmonic amplitudes defined in i|4.3l 



Table 5 

Measured variations in brightness, velocity and Ca EW in October 2008. 



mode 


frequency 


Aflux 


To,flux 


He 




T(),EW 




^vcl 




Rv 




/iHz 


% 


S 




% 


s 




km s^-*- 


s 


Mmrad"'^ 


Fl 


1210.642±0.005 


3.48±0.04 


755±1 




6.8±1.9 


150±37 


1.96±0.56 


3.6±0.6 


97±21 


13.7±2.2 


F2 


1111.806±0.005 


2.77±0.04 


310±2 




8.4±1.9 


68±33 


3.03±0.70 


2.6±0.6 


451±32 


13.6±3.0 


F3 


1224.716±0.011 


1.44±0.04 


664±3 
















F4 


1238.990±0.021 


0.66±0.04 


279±7 
















F5 


1125.730±0.023 


0.68±0.04 


151±8 
















F6 


1098.795±0.015 


1.04±0.04 


217±5 
















F7 


1292.522±0.046 


0.34±0.04 


129±13 
















F8 


1522.618±0.016 


0.94±0.04 


461±4 
















F9 


2502.286±0.014 


1.10±0.04 


249±2 
















FIO 


1999.676±0.015 


1.01±0.04 


116±3 
















Fll 


1011.818±0.021 


0.75±0.04 


294±8 
















=2F1 


2421.284±0.000 


0.30±0.04 


13±8 


2.5±0.3 














=2F2 


2223.613±0.000 


0.97±0.04 


24d=3 


12.7±0.5 














=F1+F2 


2322.448±0.000 


0.55±0.04 


223±5 


2.9±0.1 














=F1+F3 


2435.358±0.000 


0.20±0.04 


373±12 


2.0±0.2 














=F2+F3 


2336.523±0.000 


0.62±0.04 


168±4 


7.9±0.3 














=F2+F4 


2350.797±0.000 


0.20±0.04 


413±13 


5.5±0.5 














=2F4 


2477.981±0.000 


0.13±0.04 


335±19 


29.9±9.0 














=F3+F4 


2463.707±0.000 


0.21±0.04 


119±12 


10.9±1.1 














=F1+F6 


2309.437±0.000 


0.35±0.04 


200±7 


4.8±0.3 















Note: Variations are measured using the form Asin(27r/(t + To)). All Tq values are measured from BJED=2454700. 



Table 6 

Measured variations in brightness, velocity and Ca EW in November 2008. 



mode 


frequency 


Aflux 


To,flux 


Rc 


Aew 


To.EW 


ICa 


Avel 


To,vel 






/iHz 


% 


s 




% 


s 




kms~^ 


S 


Mm rad~^ 


Fl 


1210.710±0.006 


3.73±0.03 


21±1 




5.5±1.6 


334±38 


1.47±0.43 


3.2±0.7 


127±27 


11.3±2.3 


F2 


1112.606±0.008 


2.46±0.03 


192±2 




6.6±1.6 


757±34 


2.68±0.66 


3.7±0.7 


381±25 


21.6±3.8 


F3 


1224.899±0.022 


0.93±0.03 


73±4 
















F4 


1237.013±0.148 


0.14±0.03 


207±29 
















F5 


1126.067±0.035 


0.59±0.03 


25±8 
















F6 


1097.029±0.044 


0.47±0.03 


556±10 
















F7 


1293.292±0.105 


0.20±0.03 


262±20 
















F8 


1522.459±0.035 


0.59±0.03 


482±6 
















F9 


2502.317±0.018 


1.15±0.03 


322±2 
















FIO 


1999.583±0.021 


1.00±0.03 


400±3 
















Fll 


1014.329±0.020 


1.04±0.03 


605±5 
















=2F1 


2421.420±0.032 


0.65±0.03 


80±3 


4.7±0.2 














=2F2 


2225.212±0.019 


1.08±0.03 


388±2 


17.8±0.6 














=F1+F2 


2323.316±0.028 


0.73±0.03 


188±3 


4.0±0.1 














=F1+F3 


2435.609±0.194 


0.11±0.03 


166±19 


1.5±0.2 














=F2+F4 


2349.619±0.079 


0.26±0.03 


285±8 


38.4±4.9 














=2F4 


2474.026±0.091 


0.23±0.03 


167±9 


1166±409 














=F3+F4 


2461.912±0.162 


0.13±0.03 


328±16 


49.4±8.4 














=F1+F6 


2307.739±0.092 


0.23±0.03 


417±10 


6.4±0.5 















Note: Variations are measured using the form Asin(27r/(t + To)). All Tq values are measured from BJED=2454700. 
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Table 7 

Measured Ca EW and R values for the Keck and VLT observations. 



mode 


Period (s) 


Aew (%) 




{£, m) 


Kcckl 


614.20 


9.9±1.4 


3.1±0.4 


(1, 1) 


Keck 2 


817.66 


8.4itl.4 


3.1±0.5 


(1,0) 


VLTl 


615.38 


5.5±1.7 


2.0±0.6 


(1, 1) 


VLT 2 


810.76 


8.2±1.7 


3.6±0.8 


(1, 0) 



